Topological instabilities of spherical vesicles 
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Within the framework of the Helfrich elastic theory of membranes and of differential geometry 
we study the relative stability of spherical vesicles and double bubbles. We find that not only 
temperature, but also magnetic fields can induce topological transformations between spherical 
vesicles and double bubbles and provide a phase diagram for the equilibrium shapes. 
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Self- assembled vesicles were experimentally observed 
to form a variety of shapes [1]. However most equilib- 
rium vesicles are spherical. This is not surprising, be- 
cause spheres are known to minimize not only the sur- 
face energy but also the elastic energy of vesicles 0. In 
mathematics it has been proven that spheres have the 
least area to enclose and separate a given volume. The 
analogous problem of identifying the surface with least 
area enclosing and separating more than one volume was 
first considered by the Belgian physicist J. Plateau in the 
XlXth century [3] . Surprisingly, only recently it has been 
proven that a double bubble, a figure composed of two 
spherical caps separated by a flat disk (Fig. [1^), is the 
solution to the two equal volume isoperimetric problem 
in Euclidean space [4] . A general proof for double bubble 
enclosing two different volumes separated by a non-flat 
membrane was given later [5]. 

The connected spherical bubbles often observed in soap 
are certainly a realization of this 'minimum property' of 
double bubbles since surface tension is the dominant driv- 
ing force. In this paper we explore theoretically whether 
double bubbles can represent the equilibrium shape of 
vesicles when not only the surface tension, but also other 
energy terms like elastic and magnetic energies are rele- 
vant. To the best of our knowledge, a double bubble was 
never considered as a candidate to describe the shape of 
vesicles. The presence of the membrane S2 (see Fig. [1^) 
distinguishes double bubbles from dumbbells or pear- 
shaped vesicles that were extensively studied in phase 
separated systems Two adhering vesicles may look 
similar to a double bubble, but they have different topol- 
ogy, because the membrane S2 is composed of two inter- 
acting membranes [7| . Hence, when vesicles adhere there 
is no change in topology and the Euler characteristic [s^ 
remains unchanged as for two single vesicles (x = 4), 
whereas a double bubble has X = 3 [9], that is neither 
the topology of two single vesicles nor the one of a larger 
sphere (x = 2). One may therefore refer to a double bub- 
ble as an intermediate shape in the topological transfor- 
mations between one and two spheres. Such transforma- 
tions might also be relevant to describe the first stage of 
cell division. 

To study the instabilities of spherical vesicles towards 
double bubbles we start from the well known model of 
elastic free energy proposed by Helfrich ^10|] . This model 



FIG. 1: (color online) (a) The double bubble is the unique 
surface of least area that encloses two equal volumes; it is 
composed of two identical spherical caps 5*1 and 5*3 separated 
by a single membrane S2 that is flat for equal volumes; the 
three pieces meet along a circle at 120° angle [4]. (b) Cross- 
section of a magnetically deformed double bubble with semi- 
axes a > c; 2a; is the angle between two tangent planes at the 
meeting circle. 



was successful in explaining the shapes of red blood cells 
and in predicting some new non-trivial shapes of vesi- 
cles [l|, Q. However, one cannot apply without restric- 
tions a continuum elastic theory to the shapes with dis- 
continuities, since it would yield infinite energy. The dis- 
continuities can be taken into account by introducing a 
line tension with an effective angular dependence, as was 
proposed in [ll[. Similar arguments can be applied to 
a double bubble, which has a circular rim with a singu- 
larity, where three smooth surfaces meet (see Fig. [1^). 
Therefore, to take into account the rim singularity in the 
shape of a double bubble, we propose a new phenomeno- 
logical term derived from the Gauss-Bonnet theorem, 
and added to the Helfrich elastic free energy. Although 
other forms of the rim energy might be plausible, here 
we investigate the consequences of our ansatz. Assum- 
ing that intrinsic parameters of the system entering the 
phenomenological Helfrich free energy, like elastic moduli 
and spontaneous curvature, are temperature and pres- 
sure dependent [12], we propose a phase diagram for the 
studied shapes. It is intuitively clear that shapes similar 
to a double bubble may occur whenever flat membranes 
become favorable, for instance in presence of high mag- 
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netic fields that tend to orient diamagnetic molecules. 
For this case, experimental evidence of deformation of 
self-assembled spherical vesicles to oblate spheroid was 
reported in [13, 14]. Here, we find that a magnetic field 
can also infiuence the relative stability of shapes and can 
make two spheroidal vesicles less stable than a deformed 
double bubble or one larger spheroid. The latter trans- 
formation was experimentally observed for liposomes [isf 
in high magnetic fields. 

In Fig. [1^ we show the double bubble as defined math- 
ematically. The double bubble is composed of two spher- 
ical caps Si and ^'3 separated by a disk 5*2 , meeting along 
a common circle (rim) at an angle of 120° [4]. To consider 
magnetic deformations, we study the modified geometry 
shown in Fig.[T]3. Firstly, we assume that, during the de- 
formation, the spherical caps and Ss (with = Ss) 
change to oblate spheroids with semiaxes c and a parallel 
and perpendicular to the magnetic field B respectively. 
Secondly, we allow the angle between the three surfaces 
to differ from 120° by introducing the parameter uo (see 
Fig. [1)3), related to the radius of the membrane as 



V sin co -\- cos^ co 

In this paper we extend the concept of double bubble 
from the unique area minimizing surface with c = a and 
uj = 7r/3 to the set of shapes with c G (0, a], cj G (0, 7r/2). 

The equilibrium shapes of fiuid membranes are usu- 
ally studied in terms of the Helfrich model of bending 
energy [TO, [3: 

j^ei = JdS{2k{H-Hof^kK}, (2) 

where H and K are the mean and Gaussian curvatures 
respectively, and Hq is the spontaneous mean curvature. 
The first term, proportional to the bending rigidity /c, 
describes the local deviation from equilibrium curvature 
i^o, while the second term with Gaussian rigidity k in- 
fiuences only the topology and thus is often neglected. 
For the same reason, k is also not measurable if no topo- 
logical transformations occur. Here we cannot omit this 
term, because we examine different topologies. It was 
shown theoretically and confirmed experimentally that 
— 1 < k/k < for monolayer vesicles, while for bilayer 
vesicles k can be positive [12|, [13]. The values of i^o, k 
and k depend on pressure and temperature, and there- 
fore may be considered as parameters and not as intrinsic 
properties of the system. We will consider the case with 
0.7 < HqRq < 1.3, — 1 < ^//c < 1, where Rq is the radius 
of the spherical vesicle, and k = 4 - 10~^^ J ^ 10 /cbT, 
the latter being a typical value for molecular vesicles. 
Moreover, we assume that self-assembled molecules keep 
a constant density, implying the condition of constant 
surface of the vesicles S = SttRq. 

The form of the phenomenological free energy, given by 
Eq. ([2]) , is valid only for smooth surfaces like spheroids or 
the individual surfaces 81,82 and 83. In order to intro- 
duce the energy of the rim phenomenologically, we make 



use of the Gauss-Bonnet theorem J d8 K = 27rx [1], 
where x = 3 for the double bubble [9] . By integrating ex- 
plicitly the Gaussian curvature K over the three smooth 
pieces of the double bubble, namely the spheroidal caps 
{81, 83} and the fiat membrane 5*2 {K = 0), we find the 
following expression 

3 

V/ d8K = 47t{1^cosu;), (3) 

which depends only on the angle uo (see Fig. [TJd). Since 
J d8 K = Gtt, the difference with respect to Eq. ([3|) rep- 
resents the contribution of the curvature integral over 
the rim. This integral accounts for the discontinuity 
of the shape and yields the angle deficit of the rim as 
27rx — X]i=i Is d8 K = 27r(l — 2 coscj). We then assume 
that the phenomenological free energy of a circular rim is 
proportional to the square of the angle deficit (1 — 2 coscj) 
multiplied by the rim length 2'kR^ 

^rim = - 2cOSC^)'2^i?^, (4) 

where the coefficient of proportionality 7 represents the 
line tension. The value of 7 can be estimated by consider- 
ing the formation of vesicles due to edge effects and com- 
paring the free energy of the two limiting cases, namely 
a disc with Fdisc = ^tt^R^ = 4:7tjRo and a sphere with 
^sphere = 47r2A;, yielding 7 = 2k /Rq jl^. Our conjec- 
ture for the form of the rim free energy is similar 
to the line tension with an effective angular dependence 
proposed in [ijj to justify the occurrence of slope dis- 
continuities between different membrane domains. Also 
for crystals, it has been demonstrated mathematically 
in [19], that the cusps that are often present at crystal 
surfaces can be intrinsic property of equilibrium shapes 
and not necessarily due to defects. 

We begin by considering the symmetric case, namely 
a double bubble composed of two spherical caps with 

c = a = R = Rq \J 8/ (4(1 + coscj) + sin^ cj), which yields 
R\^^Q = Rq. The result for the elastic energy and 
the rim energy is = 47r/c(l + coscc;)[(l — RHq)'^ + 
k/{2k)] + 27rksin^ij{RHo)'^ and J^rim = 47r/csincj(l - 
2cosix)YR/ Rq. The cj-dependence of these terms is 
shown in Fig. [2^. The rim energy vanishes at Cc; = 
when two vesicles just touch and at cc; = 7r/3 when 
the angle between three surfaces is 2uo = 27r/3 (see 
Fig. [1]) that coincides with the Plateau's rule for soap 
films. However, contrary to soap bubbles, the elas- 
tic energy grows with uj so that the total free energy 
^tot = + ^rim docs uot ueccssarily have a minimum 
in the vicinity of co = 7r/3, depending on the relative 
contribution from J^q\. These results show that with 
this approach we can obtain the correct limiting cases. 
Figure [2)3 presents the phase diagram for the relative 
stability of spherical vesicles and double bubbles in the 
plane of the parameters {k/k)-{HoRo). Thermal fiuc- 
tuations can be taken into account via the renormaliza- 
tion of the bending and Gaussian rigidity given by (2^ 
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FIG. 2: (color online) (a) The free energy terms J^rim, ^ei and Ttot of double bubble with c = a plotted versus the angle uo. 
We choose k — and HqRq — 1 for the elastic part for illustration, (b) Phase diagram showing regions where two spheres with 
radius i?o, a standard double bubble with c — a and uo — tt/S, and one single vesicle with radius \/2Ro have the lowest free 
energy. The parameters Ho, k and k depend on temperature and pressure ^12i] . 




FIG. 3: Scanning electron microscopy (SEM) image of 6T 
vesicles on a metallic surface found after evaporation of a 
solvent. 



/Ceff = - (3/CbT/8^) log(gLx/^min) ^nd ^eff = ^ + 

(5/cB^/127r) log((7^g^j^/(7^-^). Therefore, an increase in 
temperature might be enough to induce shape transfor- 
mation. In Fig. [3] we show a SEM image of a structure 
formed by self-assembled bolaamphiphilic sexithiophene 
(6T) molecules that usually form hollow spherical vesi- 
cles [l^. The angle between two spherical parts in this 
image suggests the possibility that this structure is a dou- 
ble bubble and not two adhering vesicles. It would be 
interesting to apply other experimental techniques, such 
as TEM, to investigate this issue. 

Now consider the effect of magnetic field, which tends 
to align diamagnetic molecules and usually leads to de- 
formation of spherical vesicles towards superspheroid like 
the one measured and explained in Here, we con- 

sider the possibility that magnetic fields could affect the 
topology as well. We minimize the total free energy 



(5) 



with respect to the deformation c/a for a given value 
of CO. The elastic energy J^q\ and the rim energy J- rim are 



given by Eqs. ([2|) and (|4j), respectively, and the magnetic 
energy, assuming that the membrane 5*2 is perpendicular 
to the direction of the magnetic field, is 



mag 



2/io 



dS nl + i^Ri 



(6) 



where Ax = X|| ~ X± is the difference of magnetic sus- 
ceptibility along the long and short axes of the molecule 
and D is the length of the molecule and the last term is 
the contribution of the flat disk with radius i?m (Eq. 
Whenever the reduction in magnetic energy ^mag due to 
the presence of the membrane compensates the elastic 
energy Tq\ and the energy of the rim ^rim, the formation 
of a double bubble is favored. 

We show first the free energy, minimized over the defor- 
mation c/a, for all possible values of in a magnetic field 
chosen as 5 = 20 T (see Fig. H^). We distinguish three 
deep minima, separated by high barriers, corresponding 
to the three equilibrium shapes. By varying the value 
of Gaussian rigidity we find as the ground state: two 
spheroids for k = — 0.5/c, a deformed double bubble with 
uj :^ 0.377r for ^ = 0, and a single spheroid for k = O.bk 
with comparable values of deformation c/a. Neverthe- 
less, the presence of three separated minima means that 
the experimentally observed shapes would depend on the 
kinetics of the system and that under certain conditions 
all three shapes could exist simultaneously. In Fig. [Ud 
we illustrate the phase behaviour of the studied shapes 
at 5 = 20 T compared to the one at 5 = T in Fig. Eb. 
As we expected, in presence of a magnetic field, the area 
in the phase diagram of deformed double bubbles and 
single spheroids increases significantly, which leads to a 
larger probability of finding these shapes for the fluctu- 
ating spherical vesicles. These results also show that the 
change in topology from sphere to double bubble can be 
induced not only by temperature but also by magnetic 
fields. 
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FIG. 4: (color online) (a) Variation of the minimum of the total free energy AJT = mina;=const ^tot(c/a) — mina;=o .Ftot(c/a) 
for magnetic field B = 20 T. The thick dots on each curve indicate the local minima separated by barriers. For some cases we 
show the corresponding equilibrium shapes. For this plot the parameters are HqRq = 1, i?o = 100 nm, D — ^ nm, Ax = 10~^. 
(b). Phase diagrams enclosing three regions with topologically different shapes: thick black lines [B — 20 T) and grey thin 
lines {B — T same as Fig. [2J3). Shapes with different topologies illustrate schematically the regions of the lowest free energy. 



In conclusion, we found that the double bubble, the 
surface of smallest area enclosing two equal volumes, also 
minimizes, under certain conditions, the free energy of 
self- assembled elastic vesicles. We have explicitly shown 
that magnetic fields can be used to alter not only the 
shape and the size, but also the topology of two dia- 
magnetic vesicles. Considering the possibility of other 
topologies is useful when analyzing experimental results 



and establishing whether self-assemblies are in equilib- 
rium or not. The calculations presented in this paper 
may be thought as a first step to consider the formation 
of foam during self-assembly, like the one observed in [2l| . 
We hope that our work will provide a motivation for new 
experiments. 

It is a pleasure to acknowledge helpful discussions with 
Michael Mueger, Philippe Nozieres and Efim Kats. 
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